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Abstract. We explore and clarify the connections between two different forms 
of the renormalisation group equations describing the quantum evolution of 
hadronic structure functions at small x. This connection is established via a 



> 

■ Langevin formulation and associated path integral solutions that highlight the 

statistical nature of the quantum evolution, pictured here as a random walk 



in the space of Wilson lines. The results confirm known approximations, form 



the basis for numerical simulations and widen the scope for further analytical 
studies. 

a: 

Ph! 1 Introduction 
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A lot of effort is presently devoted to the study of hadronic systems with large 
parton densities, as revealed for instance in high energy scattering experiments. 



Of particular interest is the regime of parton saturation [|1]-|9| expected to end the 
growth with increasing energy of the structure functions which is predicted by 



linear evolution equations [p!0| , pd|] and is experimentally observed at HERA. These 
linear equations are in fact only justified in situations where the partons in the 
hadrons form a dilute system; they single out two-particle correlators (actually 
the single particle densities) as the only important contributions. In the high 
density regime however, all n-point functions have to be treated on the same 
footing, and the ensuing evolution is non-linear. Several formalisms have been 
developed recently to cope with this non linear regime p|. f\. P^H^] ■ (See also 
Refs. IKREB] for recent reviews and more references.) 



To be specific we shall focus on deep inelastic scattering where a virtual photon 
collides with a hadronic target. The crucial variable in the discussion is the large 
rapidity separation r = ln(l/x) ~ Ins between the projectile and the target; s is 
the total energy squared and x is the typical fraction of the hadron momentum 
carried by the partons seen by the photon. Large rapidity interval implies "small 
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x"and all formalisms exploit this kinematic situation to reduce the number, and 
select the type, of degrees of freedom involved in the event. 

A key simplification comes from the recognition that the basic degrees of freedom, 
in an appropriate frame and for a suitable gauge, involve only the component 
A + (x~,x) of the gauge field of the hadron, and this enters only through path 
ordered exponentials along x~ (Wilson lines): 



U\x) = Pexp [ig / dx~A+(x~, x)t a 



(Throughout, we shall use light-cone coordinates, e.g., fc M = (k + ,k~,k), with 
= (k° ± k 3 )/^/2, and k is a vector in a plane perpendicular to the direction 
of propagation of the photon, to which we shall refer briefly as the "transverse 
plane".) Thus the total cross section in deep inelastic eA (7*^4) scattering is of 
the schematic form (r = x — y, b = x + y) H,p|,|2T 



a 



7* A 



>,Q 2 



2 dz 



d 2 r 



z,r;Q 2 )\ 2 d% (l-S T (x,y)) 



;i.2) 



which naturally separates into a photon wave function t/>, giving the probability 
that the virtual photon splits into a qq pair (a color "dipole" ) with transverse size 
r, and a matrix element summarizing the QCD interaction of the dipole with the 
hadronic target: 
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S T (x,y) = —(tr(U\x)U(y))} n 



(1.3) 



This form ensues if the energy is large enough for the scattering to be treated 
in the eikonal approximation with the qq pair not deflected from its light-like 
trajectory. 

At this point, although one has achieved a reduction of degrees of freedom, the 
calculation of the average (...) in Eq. ( |1.3| ) is still a very hard problem. Note 
that this average depends on the rapidity interval r, and rather than calculat- 
ing directly the average, one can at least try to calculate its variation with r. 
This calculation has been approached by following two seemingly distinct routes. 
In the first one, one focuses on the averages of products of Wilson lines and 
obtains differential equations for their evolution with r without having to spec- 
ify explicitly the average on the target hadron |17|, 18, . Such equations form 



an infinite hierarchy coupling together correlators of Wilson lines of increasing 
order [17] (see also Refs. [[H|, 20, 24, [2(|). Remarkably, this hierarchy is coded 
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in a single functional evolution equation for an appropriate weight function 
(see also Ref. ||23|| ). The second approach builds on an explicit model for the 
hadron wavefunction which allows for a direct and explicit calculation of the 
average §,|,|,|T|, 



15,16,21,22 



In both approaches, the evolution results from a renormalization group (RG) 
procedure in which quantum fluctuations of the gluonic fields are successively in- 



tegrated out p5| 21],|22|. Reassuringly, as shown in Ref. p2 |, both approaches yield 
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identical results for the calculation of observables. Both lead to a description of 
quantum evolution of parton densities in terms of equations for probability dis- 
tributions. A complete derivation and thorough proof of these correspondences, 
however, has not been available up to now; this is the main purpose of the current 
publication. 

As we shall see, both pictures of the evolution can be vizualized as a random walk 
in the space of Wilson lines. For each value of r, there exists some probability 
that a particular value of U is realized, i.e., there is an associated probability 
distribution Z T [U] that allows the calculation of averages such as that in Eq. ( |1.3| ). 
This distribution satisfies a Fokker-Planck equation: this is the evolution equation 
of the first approach, which we shall refer too as the "[/-representation". If we 
go deeper in the details of the random walk, one finds that the value of U at 
rapidity r + dr, call it U T+ d T , is related to U T by U T+ d T = U T exp(— ia T dr), where 
a T is a random variable — the component A + of a classical colour field — whose 
probability distribution W r [a] is explicitly calculated in the second approach 
which provides also the corresponding evolution equation. We shall refer to this 
second approach as the "a-representation" . The correspondence between the two 
approaches is easy to see from the point of view of the a-representation. Indeed, 
in this representation, the matrix U at rapidity r is constructed step by step from 
some initial Uq by (right) multiplication by a matrix of the form exp(— ia v ) when 
a v is a random variable. Since we know the probability distribution of all the 
random variables a v where rj spans the rapidity interval — this is W T [a] - - one 
can easily calculate the probability Z T [U] that the random walk terminates at U. 
Furthermore, one can relate the two evolution equations. 

Further insight is gained by considering the random steps as a Langevin process, 
and by writing the average over the random variables as path integrals. Such 
path integrals provide a direct connection between the two approaches. This 
not only highlights the physical concepts driving the RG procedure with a clear 
reference to its statistical nature, but it also allows us to give simple formulae 
for the Langevin process, which may be convenient for numerical simulations. 
Investigations along these lines have actually already been carried out and will 
be available soon [^HJ. In Ref. [P5f| , a different path integral has been proposed as 
a solution to the infinite hierarchy of coupled equations established in |17| |. The 
precise relation between the results in [^] and the results that we shall obtain 
below in this paper is still to be explored. 

In the next section we shall elaborate on the physics underlying the RG equations 
of Refs |0| and |22] and recall the explicit forms of the basic equations in the 
two representaions introduced above. The formal relation between the two for- 
malisms will be outlined. In particular the characteristics of the random process 
responsible for the evolution with rapidity will put forward. Then, in section ||] 
we show that much of the structure of the evolution equations can be identified in 
the simpler setting of the random walk in ordinary Euclidean space. The last sec- 
tion considers finally the the relevant case of random walks on a group manifold. 
Conclusions are summarized at the end. 
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2 Nonlinear quantum evolution in two represen- 
tations 



As discussed in the previous section, the calculation of scattering observables at 
high energy reduces to that of averages of products of Wilson lines. This can be 
implemented in several ways, and we shall review in this section two seemingly 
independent approaches which turn out to be equivalent. Our goal here is to 
recall the main physical ingredients. The rest of the paper will be devoted to a 
deeper understanding of the relation between the two approaches by analyzing 
their common mathematical grounds. 



2.1 Wilson line operators and the U— representation 



Consider the evolution of the cross-section (|L2 ) when the rapidity gap r increases 
by an amount dr. We choose the frame so that r H rs r 3> r dipo i e , where r H and 
Tdipoie are the rapidities of the hadron and the dipole, respectively. If initially 



ct s Tdi P oie <C 1, the dipole is just a quark- ant iquark pair PB[ . When Tdi po ie is 
increased by an amount dr, so is the phase space for the emission of a gluon. 
If such a gluon is indeed emitted, then the interaction of the dipole with the 
target will involve the independent scattering of the quark, the antiquark, and 
the emitted gluon, off the colour field A(x~,x) of the target. This is the source 
for the change dS T oc a s dr in the S'-matrix element S T = 4-(tr(C^ Z7 y )) T , for 
which the following evolution equation is obtained: 

d T {tr(UtU y )) T = 

a s f j2_ (x-y) 



2tt 2 



/ [x-z)\y-zf ( N MUSP V ) - tr(UlU z )tT(UtU y )) T , (2.1) 



where all the Wilson lines are in the fundamental representation, and U x = U (x) . 
Eq. Q2.1p has been originally derived by Balitsky JT7 . 



While this equation for S T (x,y) involves only averages of products of Wilson 
lines, this is not a closed equation: the 2-point function is related to a 4-point 
function. However a closed equation is obtained when one allows N c to take 
arbitrary large values. Indeed, in this large N c limit, the product of traces in its 
r.h.s. factorizes, 

(tr(UiU z )tr(UiU y )) T — <tr(^,)> T (ix(UlU v )) r for N c - oo, (2.2) 

so that eq. (|2J]) reduces indeed to a closed equation for S T (x,y). In this form, 
the equation has been independently obtained by Kovchegov |18|], within the 



Mueller's dipole model It has attracted much interest recently [|30|-|33||.(See 



also Ref. |20| for still a different derivation, which is summing "fan" diagrams.) 

In general, eq. ( |2.1| ) is just the first step in an infinite hierarchy of coupled evo- 
lution equations for the correlation functions of Wilson lines. In principle, these 
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equations can be all obtained in the formalism of Ref. [17], but, as recognized in 
Ref. [H]], they are equivalent to a single evolution equation for an appropriate 
generating functional Z T [U] giving the probability that a certain field U x is real- 
ized in a collision with rapidity gap r. In particular the averages of products of 
Wilson lines are obtained as: 



(U£ • • • L© T = JmU)] US ■ ■ ■ U$> Z T [U] . (2.3) 
Here, t/ (t) is a generic notation for either U or W (in any representation of 



SU(iV)), and d/j,(U) denotes the group invariant measure [5B, 57J (see also Ap- 
pendix). Eq. (|2.3| ) is consistent with Balitsky's equations if Z T [U] obeys the 
following equation: 

drZ T [U] = ^lxty[U)V b y Z T [U), (2.4) 

where we employ a summation and integration convention for repeated indices 
and transverse coordinates - - in the reminder of the text we will only write 
integral signs where confusion might arise. is a Lie derivative which generates 
translations on the group manifold (see Appendix), 



,ab 



XZP] = \J^i - Vl U »)0- ~ UlU y )r (2-5) 



with U in the adjoint representation, and 

= ( x -z)-(y- z) 
xyz ~ {x-zf{z-yf 1 } 

X is of the form e^e and hence positive definite. More precisely 

xtiU] = I d 2 ze ac ' l (x,z)e bc ' l (y,z), (2.7) 

which, according to our conventions, we could have written shortly as Xxyt^l = 
e< xz eb yz- The "square root" factor in Eq. ( |2.7| ) is given by 



1 (x — z 



e 



It is a matrix in color indices and a vector in transverse coordinates. With the 
help of 

<[U] = l^UtPl ( 2 - 9 ) 
Eq. ( p.4| ) takes the familiar form of a Fokker-Planck equation 

drZ T [U] = \ V a x V b y {xt v [U\ Z T [U}) - VI {a a x [U\ Z T \U}) , (2.10) 
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albeit a functional one, involving fields with values on a curved manifold. As 
for all Fokker-Planck equations, there is an underlying random process whose 
origin lies in the renormalisation group operation itself, and the elimination of 
degrees of freedom which accompanies a change in r. In fact, this aspect will 
be made the centerpiece of the formulation by recoding Eq. ( |2.4j ) as a random 
walk governed by a Langevin equation. However this short presentation does not 
make the elements of the random walks immediately visible; it does not either 
exhibit the important role of the background field of the target in determining 
the properties of the random walk. Those aspects are better seen with a picture 
of the target wave function which utilizes the field a rather than U as the basic 
random variable, a description to which we will turn now. 



2.2 The ^-representation 

The model of hadron wavefunction that we shall use rests on the effective theory 
developed in Refs. 0, |], [15|, ^2] , which we briefly summarize (see Ref . |25j for 
details) 

Because of the large rapidity gap between the dipole and the hadron, the latter 
appears to the dipole as a Lorentz contracted colour source for the small-x gluons 
to which it couples. We are then led to separate the hadron constituents into hard 
ones with large longitudinal momenta (k + ^ xP + ), and soft ones with momenta 
k + ^ xP + (with P + = the total longitudinal momentum of the hadron). The hard 
partons include the valence quarks, as well as the partons created in the quantum 
evolution down to longitudinal momenta xP + and whose density increases with 
1/x. Together, these hard partons constitute a high density system of colour 
charges which, during the interaction with the dipole, can be regarded as frozen 
and constitue the source p a (x~ , x) of a classical field representing the soft gluons. 

This classical field is assumed to be the solution of the Yang-Mills equations: 
(D u F Ufl ) a (x) = p a (x~ , x) , where the source p a (x~,x) is time- independent, 
i.e. independent of x + . This assumption follows from the observation that the 
internal dynamics of the hadron is frozen during the relatively short interaction 
time Ax + ~ 2xP + /Q 2 with the dipole. The classical field has just one inde- 
pendent component, which in the covariant gauge d^A^ = is the component 
x) = a a (x~, x), with -d\u a {x~, x) = p a (x~, x) . 

The source p a , or equivalently the classical field a a , is treated as a random vari- 
able, with some gauge- invariant probability density W r [aj , which depends upon 
t = ln(l/x) |5],[IJ|. In order to compute some observable 0[a] (to leading log 



accuracy), one first evaluates it with the classical solution, and then average the 
result with the weight function W T [a] : 

(0) T = J [da]0[a]W T [a]. (2.11) 

This peculiar averaging, somewhat reminiscent of the mathematical description 
of amorphous materials like spin glasses, has been dubbed Colour Glass Con- 
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densate [22]. What is being done here is analogous to the Born-Oppenheimer 
approximation: the original quantum average is replaced by a classical average 
over the various possible configurations of the "frozen" degrees of freedom. 

The probability distribution W r [a] is not known directly, but its variation corre- 
sponding to freezing partons in the rapidity window (t, r + dr) can be computed. 
These partons have longitudinal momenta bA < k + < A, where A is the mo- 
mentum scale corresponding to r (i.e., r = ln(P + /A)) and dr = ln(l/6). The 
elimination of these degrees of freedom induces new correlations at the softer 
scale bA, which, to the accuracy of interest, can be accounted for by shifting the 
original field a a (x~,x) by a random quantity a^(x) dr with only 1-point and 
2-point correlation functions: 

(a a T (x)) A , {a a T (x)a b T (y)) A , (2.12) 

where the subscript A is here to remind us that the averages depend on the 
background field a(x~,x) at the scale A. 

At this place we need to digress on a subtle aspect of the problem which has to 
do with the longitudinal structure of the field. The random field a^(x) is a two 
dimensional field which is in fact obtained by integrating over x~ a fluctuation 
5a®(x~ , x): 

a a T (x)dr= / dx-5o£{x-,x). (2.13) 



It is shown in |£2| that the x dependence is of the form: 

g— ibhx~ g— iAx - 

6a°(x~,x) oc 0(aT) = , (2.14) 



x~ 



so that the integral over x in Eq. (|2.13|) produces a factor ln(l/6) = dr. 
Eq. ( |2.14j) shows that the shift in the classical field resulting from integrating 



out gluons in the chosen layer of k + has support at 1/A ^ x~ ^ 1/&A, and there- 
fore essentially no overlap with the original field a a (x~ , x) at the scale A (whose 
support is x~ iS 1/A). This is an important feature which allows us to reconstruct 
step by step the Wilson lines involved in the calculation of high energy scattering 
matrix elements. Thus, for instance: 

Ul +dT (x) = e i9fdx-Sa«(x-,*) [ft (a;) = e&«*ra?(«0 jji^ ? ( 2 .15) 

since all the x~ in 8a® (x~, x) are larger than those contained in the Wilson line 
Ul(x). This being recognized, it is not essential to keep trace of the detailed x~ 
dependance of the field exhibited in eq. ( 2.14 ), and we shall work with integrated 



quantities such as a T (x): the ordering in r takes care of the initial ordering in 
x~ . We are now in position to introduce the equation obeyed by the probability 



distribution W T [a] |22| 



™ = \ j^m (xiWWVH) - ^ KWWM) , (2-16) 
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where X^tM an d &x[ a \ are gi ven by Eqs. ( j2.5|) and ( |2.9| ) respectively, in which 
U and W represent Wilson lines in the adjoint representation, built from the 
classical field using Eq. (|2.15| ) iteratively from tq to r. In Eq. ( |2.16| ), the functional 
derivatives are taken with respect to the random variable of the last rapidity bin, 
i.e. a T , which, according to the previous discussion, is the only one involved in 
going from PV r [a] to W r+ d T [a]. 



2.3 Path integral solutions and Langevin representations 
of the RG 

As emerges from the discussion in the previous subsection, the two evolution 
equations in the U and the a representations, summarize a random process by 
which the Wilson lines are built. This process is more directly visible in the a 
representation: there we see clearly that at each step in rapidity, the Wison line 
is multiplied by a random matrix close to unity. The Wilson line at rapidity r is 
given by 

U}(x) = Texp^ig^dria^x^Ulix), (2.17) 

with, as already mentioned, the inital ordering in x~ (cf. eq. ( p. . 1|) ) replaced by an 
ordering in r. With this explicit representation at hand, one has already indica- 
tions about the relation between the two approaches. Indeed, the Lie derivative 
which enters the equation for Z[U], and which produces a "translation" of U, is 
equivalent when acting on U such as ( |2.17| ) to: 



which is precisely the variation generated in the a representation. In fact one 
can write an explicit formula relating the two probability distributions Z T [U] and 
W T [a\: 

Z T [U] = J [da] 8 G (U\ TeH T *)^) W T [a] . (2.19) 

This also immediately allows us to translate the equations for the two into each 
other, giving complete equivalence of the approaches. 

Thus the basic problem of the evolution at small x is that of a random walk in 
a space of two dimensional fields with values in SU(N). The group nature of the 
variable introduces several complication in the random walk which will be dealt 
with in section 4. But the general structure of the problem can be uncovered 
without dealing with these complications. This will be done in the next section. 
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3 Preliminaries: Brownian motion in flat space 



As recalled in the previous sections, the evolution with rapidity of scattering 
observables can be described by two distinct functional equations. Both can 
be interpreted as diffusion equations, but they have rather different structures 
and control the evolution of quite different object: a probability distribution on 
group valued fields in one case, a functional of gauge fields in the other. It is the 
purpose of this section to clarify the deep relation between these two formalisms 
by showing that they have counterparts in the description of ordinary Brownian 
motion. 

We shall consider the Brownian motion of a particle in a D-dimensional Euclidean 
space. As we shall see, the probability to find the particle at point x at time r, 
which we call P(x, r), is related to Z T [U], with the coordinate x corresponding to 
U. Similarly one may introduce a probability functional Wt[v] over the random 
paths, where the paths are parametrized by the velocity v of the particle at each 
time step. This functional W T [/u] is the analog of W T [a]. It is related to P(x,r) 
by an equation analoguous to Eq. ( |2.19| ): 



P(x,r) = J[dv]S {D) (x-x[v])W T [v] (3.1) 



with 



x[v] = x + / dr'v(r') 
Jo 



(3.2) 



with T an arbitary time greater than r. Note that in the r.h.s. of Eq. ( |3.1| ) 
all the time-dependence is carried by the weight function W r [v]. Of course, the 
paths which are physically relevant for computing P(x, r) terminate at r: this 
information is encoded in the weight function W r [v]. 

The probablility distribution P(x,t) obeys a Fokker-Planck equation which, in 
analogy with Eq. ( |2.19| ) satisfied by Z T [U], we choose of the form 

P(x,t) = -HP(x,t) (3.3) 

where the dot denotes the time derivative, H is the following "hamiltonian" : 

H = -d a d^ X a(3 {x) -cx a {x) , (3.4) 

(with d a = d/dx a ) and x al3 { x ) is of the form: 

X a/3 {x) = e a ^{x)e^{x). (3.5) 

Thus x a ^{ x ) is symmetric and positive semi-definite, and therefore invertible. 
We shall also assume the following relation, analogous to Eq. (fX 



a a (x) = -d^ X al3 (x) (3.6) 
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which, in particular, guarantees that H in Eq. 
satisfied by W T will be found to be: 



is hermitian. The equation 



d T W T [v) 



5 

5v<i 



1 5 



X^(x)W T [v]) - a a (x)W T [v] 



(3.7) 



Note that the Fokker-Planck equation fl3.3|) -( j3!4|) for P(x, r) follows indeed from 
Eq. ( |3.7| ) for W r [v] together with the relation ( |3.1| ), as can be checked by per- 
forming some integrations by parts with respect to v in Eq. (|3.1|) and then using 
the identity: 

— 5{x - x[v\) = -d a 5(x - x[v\). (3.8) 

Underlying all these equations is an elementary random process which controls 
the random walk. This elementary process is best seen at the level of a Langevin 
equation, from which in fact all equations can be derived. We shall make this 
Langevin equation the starting point of our discussion. 



3.1 Langevin equation and path integral 

Consider then a particle undergoing a random walk in a D-dimensional Euclidean 
space. Its coordinates x a (a, (3 = 1, • • ■ , D) obey the following Langevin equation: 

x a {r) = a a (x) + e al3 (x)u l3 (T), (3.9) 

where the dot denotes time derivative, and a a (x) and e af3 (x) are given regular 
functions of the coordinates characterizing the medium in which the particle 
propagates; finally, v a {r) is a Gaussian white noise (see below). We shall also 
use the more compact notation: 

x(t) = a(x) + e(x)u(r), (3.10) 

where x, er and v are D- dimensional vectors and e is a D x D matrix. 

Strictly speaking, Eq. ( |3.9| ) is formal (the trajectory x(t) of the random walk 
is not differentiable) and gets meaningful only with a discretization prescription. 
Here, we shall focus on the following discretization: 

x i x i-l /o i 1 \ 

= <Tj_x + ej_ x i/j, (3.11) 



and only briefly mention later other possible choices (see also Refs. |34|,|35|]). In 
Eq. ( |3.11| ), e is the length of the time step, x,i = x{ri) with n — ie, i — , 1 , • • -n, 
and ceo = x(to) is a given initial condition. We have also introduced the simplified 
notation <Tj = a(xi), ej = e(xi), to be used throughout. 
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The Vi 's are Gaussian random variables with zero expectation value and local 
2-point correlations: 



(z/Vf) = -5 a %. (3.12) 
J e 

In other words, the probability law for the noise variables Vi is a normalized 
Gaussian distribution: 

dP{Ui) = (^) D/2 e-^d D Ui, (3.13) 

the same for each time step, and for n time steps it is the product of n such 
distributions. 

The particular discretization ( |3.11| ) of the Langevin equation is such that, for 
fixed random variable linearly related to i^. This makes it easy to 

calculate the probability P e (xi\xi-i) that the random walk brings the particle 
from position ccj_i to Xj in one time step, and this for any value of e. (We shall 
often refer to quantities such as P e (xi\xi~i) as "elementary propagators".) We 
have: 

P e (cc i |cc i _ 1 ) = J dV{h>i) 5 D (x,i - Xi_ x - e(<7i_i + e^i^)). (3.14) 
The integral over Vi is easily performed, and yields: 

P^Xilxi-i) = ^ et ^~)D/2 2 e-^-^- 1 '^- 1 )^^-^- 1 -^-^ , (3.15) 



where x 1S obtained from e via Eq. (|3.5| ). Thus, for fixed Xi-i, the probability to 
find the particle at position cc, is a gaussian centered around Xi = x^\ — eo"i_i, 
with a width eXi-i- 

The coordinate a; n of the particle after n time step is also a random variable. 
The average of any function of x n , say f(x n ), can be obtained by first solving 
Eq. (|3.11| ) for a given realization {ui,u 2 , - ■ ■ ,v n } of the random variables, and 
then averaging over all such realizations. Equivalently, one may first determine 
the probability density P{x,r) to find the particle at point x at time r = ne 
knowing that it was at Xq at time tq = 0. The average of f{x) may then be 
obtained as 



(f(x n )) = (f(x)) T = J d D x f(x) P(x, r) (3.16) 
The probability P(x, r) is given by: 

/n 
Y[dV(ui) 8 D (x-x n [u]), (3.17) 
i=i 
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where x n [u] is the solution of Eq. ( p.llj) for n time steps, and depends therefore 
on all the i/j's with i < n. Using repeatedly Eq. ( 3.14Q one can easily show that 
the integral above can be rewritten as: 

/n— 1 n— 1 

II d ° X i P&\*n-l) II Pe(*i\*i-l)- ( 3 - 18 ) 
1=1 1=1 



By inserting in this expression the formula ( 3.15 ) for the elementary propagator, 
we finally get: 



with 



= n ji^m ri( det ^ 



-1/2 



1 (2.^ 11 



(3.19) 



(3.20) 



and 



■4W-jE( 

2=1 



(Ti-l )Xi-l 



(3.21) 



The average (f(x)) T is then easily obtained: 

(/(aj)) T = Jd D xf(x)P(x,r) = J[dx] f(x n ) e~ A[x] 



(3.22) 



An alternative expression for P{x,t) is obtained by integrating over auxiliary 
momentum variables which allows one to get rid of the determinant detXi-i in 
the measure and to reexpress A in terms of the matrix Xi-i rather than its inverse: 



P(x, 



with 



n 



'jD d °Pj 

ax, 



3 (0^\D 



(2tt 



5(x - x n ) exp < - y^Pj ■ (xj - Xj-i) + eCj] 

3=1 

(3.23) 



(3.24) 



The formulae (|3.19|) and (|3.23|) express the probability P(x, r) as a sum over paths 
specified by the set of coordinates {x , Xi, ■ ■ ■ £C n _i, x n }, where x is a given ini- 
tial condition while all the other coordinates are random variables. A noteworthy 
feature of this expression is that, at step i, the functions x( x ) an d have to 
be evaluated at the point X{-i attained in the preceding step. This property is 
of course intimately related to the special discretization of the Langevin equa- 
tion, Eq. ( p. 11| ). One could naively expect that, in the limit e — > 0, £Cj_i can be 
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replaced by x^ in x 01 °- This is not so however, and it is important to keep 
the arguments of \ and a as written in Eqs. ( |3.19j ) or ( |3.23| ). This is because 
the typical paths contributing to the integral are such that \xi — Xi-\\ ~ y/e. 
Thus, the replacement of x^\ by Xi within \ or a modifies the exponent in the 
integrand of these formulae by terms of the same order of magnitude as the terms 
that have been kept. Therefore, writing Eq. ([3.23 ) in continuum notations, i.e., 
as 



P(x,t) 



J [ da3 ] j [ d P] ex P |^ [ - iP ' (* - <7( x )) ~\p- x(x) ■ p] 



(3.25) 



is ambiguous. 



3.2 The Fokker- Planck equation 

The path integral in Eq. (|3.19|) summarizes in a compact formula the averages 
performed at each step of the Langevin process. In the limit e — > 0, it is possible 
to obtain for P(x, r) a differential equation which involves only the local prop- 
erties of the medium in which the particle propagate, i.e. the functions \ and a, 
and with no explicit reference to the specific elementary random process. This 
equation is a Fokker-Planck equation. 

To obtain this equation we note first that Eq. ( p. 18 ) allow us to write: 



P(x,t) = J d x„_iP e (aj|aj n _i)P(aj n _i,r n -i). (3.26) 

In order to calculate the integral over ac n _i, we use the fact that the propagator 
P e (x\x n _i), considered as a function of x n _i, is a gaussian which, when e — > 
0, is strongly peaked near x (see Eq. ( |3.15| )); it follows that, when e is small, 
one can get a good approximation of the integrand of Eq. ( |3.26| ) by expanding 
P(x n -i,T n -i) in the vicinity of x. To do that systematically, it is in fact easier 
to go back to Eq. ( 3.14j) , and expand the 5-constraint, keeping all terms that will 



contribute to 0(e) after averaging over the noise: 

5 D (x - x n ^ - eK_a + e n ^v n )) = +°»-™> 8 6 D (x - x n ^) 

^{l-{(T n _ l + e n _ l u)-d + \{{a n _ l + e n ^u)-d) 2 y D {x-x n _ 1 ) , (3.27) 

where d a = d/dx a . Averaging over v n leaves us with the expression 

P e (aj|aj n _0 = (l - ea a n _ x d a + | xl\d a d^ 6 D (x - x n _ x ) + O^' 2 ) . (3.28) 

The first term obtained after plugging this formula into Eq. ( |3.26| ) is simply: 

P{x, T n -x) = P{x, r - e) w P{x, t) - eP{x, r), (3.29) 
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where P(x, r) = d T P(x, r). In the remaining terms, we use identities like 



a a {x n _ 1 )d a 5 D {x - £c„_x) P(cc n _i, r n _i) = d a (S D (x - x n ^ 1 )cr a (x)P(x, r n _i)), 

(3.30) 

to deduce that, in the limit e — > 0, P(a;, r) obeys the following equation: 



P(«, r) = <9 a [fl" (-x^(*)P(«, r)J - a a (x)P(x, r) 



(3.31) 



3.3 Alternative discretizations 

At this point, one could reverse the logics and consider the Fokker-Planck equa- 
tion as our starting point. Then, path integral representations for the solution 
P(x, t) could be constructed in a standard way, for instance by exploiting the 
formal analogy between Eq. ( 3.31| ) and the Schrodinger equation in imaginary 



time. Because the coefficients x( x ) an d <?(x) i n this equation are coordinate de- 
pendent, different orderings of the derivatives within the Hamiltonian ( |3.4|) will 
lead to different path-integral representations (possibly involving derivatives of x 
and a), which become however equivalent in the limit e — > 0. 

To illustrate this, let us consider three alternative expressions of the Hamiltonian 
( |3.4j ) which are obtained by successively commuting the derivatives with x-, an d 
using Eq. (|3.6| ) to simplify the results: 

H = - (i d a d f3 X al3 (x) - d a a a {x)) (3.32a) 
= - (\ x al3 {x)d a d /3 + a a (x)d a ) (3.32b) 

= ~\ {d a , X aP {x)}} + I {F&^ix)). (3.32c) 

The form in Eq. ( |3.32a| ) is what naturally came out of our original Langevin 
discretization, the other two represent different but clearly equivalent order- 
ings of momentum and coordinate operators; the term involving a double anti- 
commutator in Eq. (|3.32cj) is in the so-called Weyl-ordered form. By using the last 
two expressions, Eqs. ( |3.32b| ) and ( 3.32c ), within a canonical construction of the 



path-integral, one obtains expressions similar to Eq. ( |3.23|) , but with a and x eval- 
uated at Xi (for Eq. (|3.32b| )) and, respectively, at the midpoint Xi = (xi+Xi-i)/2 
(for Eq. ( |3.32c| )). The only difference with respect to Eq. ( |3.23| ) consists in the 
replacement of in that equation by the alternative expressions 

d = ~ X aP (^)ptpf + ip> a {xi) (3.33a) 



■X 



aPf Xi)p^ + ^d^x^ixi). (3.33b) 



It can be verified that the expressions of P(x, r) obtained from Eqs. ( p.33[ ) satisfy 
indeed the Fokker-Planck equation ( |3.31|) in the continuum limit. 
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Taking this one step further one can translate these differently discretized path 
integrals into differently discretized Langevin systems. Let us illustrate this for 
the case of Eq. ( |3.32b|) with (|3.33a| ) as its discretization prescription for the path 
integral. One can verify that its elementary propagator can be written as the 
following integral over an auxiliary random noise 

P t {xi\xi-i) = J dV(b>i) 5 D (x,i - aji_i - e(-(Ti + eiVi)). (3.34) 

The 5 function in this equation imposes the proper discretization of the corre- 
sponding Langevin equation. Note that this has cr and e evaluated at X{ instead 
of Xi_i and the sign of the cr term is changed compared to Eq. (ft.ll ). For the 



Weyl ordered form, the cr term would be completely absent but additional terms 
would arise from eliminating the contribution e-^l 9 " 913 *™ 13 from the measure in 
order to recover the structure of Eq. ( |3.34| ) necessary to read off the Langevin 
equation. 

From here it is evident that continuous notations, which are not able to keep track 
of details of coordinate dependence, are ambiguous, both for path-integrals and 
Langevin descriptions. The naive continuum notation for the Langevin equation 
corresponding to Eq. ( |3.34| ) would leave us with the sign of the cr term reversed 



compared to that in Eq. (|3.14j) , while "forgetting" the mid point prescription of 
the Weyl ordered version would lead to yet another form. The same holds for 
path-integral expressions involving the differently ordered C given above. 

The one common denominator of all these equivalent, but differently discretized, 
path-integrals and Langevin systems is the fact that they correspond to one com- 
mon Fokker-Planck equation which is perfectly meaningful in continuum notation. 



3.4 The velocity representation 

Until now we have focused on the x n as the important random variable. There 
are situations however where one needs to calculate observables which depend on 
the entire path, and not on the end point alone. Thus, we are led to consider 
distribution functions over the paths. 

To proceed with the construction of such a distribution, it is useful to perform a 
change of variables, and to better specify the paths. 

Let us then consider as new random variables the velocities: 

Vi = — — = <Ti-i + ei-\Vi , (3.35) 
e 

for i = 1, 2, • ■ ■ , N, where T = Ne is a fixed time much larger than all times of 
interest. The coordinates of a path are easily expressed in terms of these: 

i-i 

Xi_ x = x + e Vj = Xi^lv]. (3.36) 
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Note that this change of variables is not a linear transformation since, as is evident 
from Eq. ( |3.36[ ), the coefficients <Tj_i = (r(xi-i) and ej_i = e(xi-i) in Eq. ( p.35| ) 
for Vi involve the velocities Vj at all the steps j = 1,2, — 1. However, 
<7j_i and ej_i do not depend upon itself, so that the corresponding Jacobien 
is rather simple, being the determinant of a triangular matrix. Accordingly, the 
probability distribution for the v^s is immediately obtained from that of the cCj's 
in Eq. (|3TT9| ): 

dP[v] = [dv]e~ A[v] , (3.37) 

with 

^ = II (f) dL>v * ( det X<-i) " 1/2 , (3-38) 

i=l 

and 

v 



- fi-i) («« - ffi-i) • (3-39) 



We can write dV[v] in the convenient factorized form: 

A? 



dP[«] = JJd%Pi[u]. (3.40) 



i=i 



Although this factorization property is a major advantage of the velocity rep- 
resentation, it should be kept in mind that the notation is deceiving since the 
different factors do not correspond to independent probabilities: the probabil- 
ity density Vi[v] for the velocity at step i depends also upon the velocities V\, 
V2, ■ ■ 'Vi-i. Giving this, a convenient way to check the normalization condition 
j dV[v] = 1 (or to do any other similar calculation) is to perform the integra- 
tions over the v^s in decreasing order of i, from i = N down to i = 1 (notice that 
Eq. ( ^.37| ) is truly a Gaussian in v N ). 

Let us turn now to the specification of the paths. The paths which contribute to 
P(x, t) have a number n = r/e of time steps, and the probability w n (vi, ■ ■ • , v n ) 
that a given such path is realized in the random walk is given by Eq. ( 3.40|) in 
which iV is replaced by n, that is 



W n (Vi, 



v n ) = Y[d D v i V i [v). (3.41) 



Now we want to be able to compare such distributions for different numbers of 
time steps, say n and n + 1, and eventually obtain a differential equation in 
the limit of infinitesimal time steps. Clearly one cannot directly compare w n and 
w n+ \ since these have different numbers of variables. To allow such a comparison, 
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one possibility is to consider function of n + 1 variables in which the last 

one, v n+ i, is set equal to zero. A normalized distribution which achieves this goal 
is w n (vi,--- , v n )5(v n+ i). More generally, we shall consider paths with a fixed 
number N of time steps such that the velocities associated with the last N — n 
steps vanish, while the velocities of the first n steps are random variables chosen 
with the distribution given above. That is, we define the following probability 
distribution: 

n N 

wM = n vm n * d m- ( 3 - 42 ) 

i=l j=n+l 

This distribution depends on n, that is on time, and satisfies an evolution equation 
that we now determine. To this aim, we use Eq. ( |3.42| ) to write: 

W n [v]-W n ^[v] = mpJ ( n 5 D ( Vj ) \ (v n [v]-5 D (v n )y (3.43) 

\i=l / \i=n+l / 

Then we note that T^w], considered as a function of ev n , is a gaussian peaked 
at v n = <J n -i, and proceed as for Eq. ( |3.28 ) to write: 

In writing such an expansion, we anticipated on the fact that the relevant variable 
is actually ev rather than v. it is in the variable ev n that the width of the gaussian 
is of order e; furthermore, as we shall see shortly, ev n is the variable which is 
needed to turn W n [v] into a functional of the continuous function v(t). 

Using the expansion Q3.44 ) in Eq. ( |3.43|) one easily obtains: 

wM-w»-M - (y^A - ^<-)^-M- (3.45, 

Recall that Xn-i and & n -i involve only the velocities Vi with i < n — 1, cf. 
Eq. (|3.36|) . However, since the function W^ifv] on the r.h.s. is proportional to 
Yli> n ^ D ( v i) one ma y re pl ace x n-i by x[v] in Xn-i and <r„_i, with 



JV 

x[v] = x + e^Vi , (3.46) 

i=l 



and consequently Xn-i by x( x ) and cr n _x by a{x). The continuum limit e — ► 0, 
— > oo, and n — > oo with ne = t and A^e = T fixed, can then be taken. In doing 
so, we use the formula 



e ov n ov{t) 
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which transforms the partial derivative w.r.t. v n into a functional derivatives 
w.r.t. v(t). We thus obtain the following functional differential evolution equa- 
tion for W T [t/l : 



d_ 



WJv] 



_5_ 

5vi 



1_5_ 



X ^(x)W T [v]) - a a (x)W T [v] 



(3.48) 



where v t = v{r) and x = x[v] given by Eq. (|3.2| ). 

Once the probability distribution over the paths W T [v] is known, it is of course 
easy to calculate probability P(x,t) according to Eq. ( |3.1| ). The only quantity 
which depends on r on the r.h.s. of this equation is W T . Taking the time derivative 
and using the equation ( |3.48| ), one easily reconstructs the Fokker-Planck equation 



( 3.31 ) satisfied by P(x,r). 

One may also rewrite the probability distribution ( |3.42| ) in the alternative dis- 
cretization in which a and \ are evaluated at Xi, as in Eq. ( |3.33a| ). As in Sect. 



3.3 , this requires using the relation (|3.6|) between \ an d cr. One thus obtains: 



WJv] 



N 

n 



2tt 



D 



d D pi exp 



N 



V; 



X exp {~ {Pi ■ X(xi)- Pi + 2i pi -aixi)^, 



(3.49) 



which can be recast in continuum notations as follows: 
W T [v] = J [dp] exp | -i jf dr) p v ■ v v | 

x exp \-\J dr l (l>v ■ X(x v ) -p v + 2ip v - a(x Ti 
with x v = x v [v] given by the continuous form of Eq. (|3.36|) : 



xJv) = x + / dr'v(r' 



(3.50) 



(3.51) 



Although, for the reasons discussed in relation with Eqs. ( |3.25| ) and ( p.34| ) above, 
such continuous notations can be ambiguous, they can be convenient. For in- 
stance, when used with care, they allow simple formal operations, like functional 
differentiations. In particular, it is easy to verify in this way that Eq. ( p.50| ) 
satisfies equation ( ^.48| ). To this aim, it is convenient to use the alternative form 
of Eq. (|3.48| ) which is obtained after commuting the two functional derivatives in 
its r.h.s. through \ (i- e -> the analog of Eq. ( 3.32b| )). In this check, the following 
identity is also useful: 



8 



dr}F{x v ) = Oir-n) / dr]d a F(x 



(3.52) 



Tl 



which vanishes for t% > t. 
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4 Brownian motion in the space of Wilson lines 



After all these preparations, we finally come to the study of the Brownian motion 
in the space of Wilson lines, which is our main interest here. The forthcoming 
discussion will in fact follow closely that in Sect. ^: as we shall see, the evolution 
equations fl2.4| ) for Z T [U] and ( |2.16| ) for W T [a], are the close generalizations of 
Eqs. ( p.31| ) and ( |3.48| ), respectively. 

The Wilson line U(x) is a field, with the coordinate x leaving in the trans- 
verse plane. In what follows, we shall often omit these transverse coordinates in 
intermediate formulae, in order to simplify the notation and focus on the true dif- 
ficulties. These coordinates will be reinserted in the final formulae, or whenever 
ambiguities may arise from omitting them. Similarly, we shall alternate between 
discrete and continuous time: most of the formulae are established for discrete 
times, and while many of them have a straightforward continuous extension, one 
should keep in mind that the use of continuous notations may hide ambiguities, 
as emphasized in the previous section. 



4.1 Langevin equation and the a- representation 

The physical random variable in the evolution is the elementary contribution 
a T {x) to the classical field of the target arising from integrating out gauge field 
fluctuations in the the rapidity strip [r, r + dr\. Such a contribution changes the 
Wilson line according to (r = ne, a T dr = a n e): 

Ui = {/i-ie"^, U} = e^ ta Ul v (4.1) 

A path in the space of U fields is thus defined as the set of values {U\, • • • , Uq}, 
with : 

C/t = e^Ul^ = ■■■ = e^V 6 "™- 1 • • • e ieai U j , (4.2) 

In continuous notations, and with spatial coordinates included, a path is given by 
Eq. ( |2.17| ) in Sect. [2.3| . Note that in the previous formulae, and in most formulae 



to come in this section, we set g = 1 for simplicity. To reintroduce the coupling 
constant g in the final results, it is enough to reinsert it in the exponent of the 
Wilson lines. 

The random variables ot%(x) at step i are selected according to a Gaussian prob- 
ability distribution with the following characteristics: 

«(*)> = aU{x) , «(*) a\{y)) = \ xf^x, y) , (4.3) 

These quantities are calculated in terms of the Wilson lines obtained at the step 
2 — 1, according to Eqs. ( |2.5|) and ( |2.9|) in Sect. 0. 
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Equivalently, the variables a^(x) may be obtained by solving the Langevin equa- 
tion (with given by Eq. (|2.8|) with U = Z7j_i): 



= ot_i{x) + J z ett(x, z) v}\z) , (4.4) 
where uf(x) is an auxiliary Gaussian white noise: 

{vf{x)) = 0, (vf{x)vf{y)) = - S^^S^ix - y), (4.5) 



By averaging over the noise the solution of Eq. ( |4.4j) one indeed recovers the ex- 
pectation value and correlator of the random variables oif(x) as given by Eq. (|4.3|) . 

The Langevin equation ( |4.4| ) plays the same role for the random walk on the group 
manifold as Eq. ( |3.35| ) for the random walk in a flat space, with the gauge field 



a%(x) playing here the role of the velocity Vi. But in contrast to what happens 
in the simple example of Sec. pT4] , the relation between the "coordinate" U{ and 
the "velocity" a, is here non trivial, because of the non-trivial geometry of the 
group manifold. However, as long as we characterize the random walk by small 
displacements in the tangent space, the special geometry of the group manifold 
plays essentially no role (it only enters through the dependence of a and \ on the 
Wilson lines) . The mathematical analogy with ordinary Brownian motion in the 
velocity representation is then complete, and the formulae obtained in Sect. |3]4 
have an immediate translation to the present problem. 



Thus the weight function W n [a] is given by the analog of Eqs. (|3.42|) and (|3.40| ), 
namely: 

n N 

W n [a]=JlVi[a] J] 5 N> ( 4 - 6 ) 

i=l j=n+l 

where 

Vi[a) = (^) D/2 (det X .-i)" 1/2 e- A[a] , (4.7) 



AN = / (a, - a^iyjx£&{<*i - *i-i) h y , (4.8) 

Jx,y 

and 

8[<*i] = nn *(«?(*))■ ( 4 - 9 ) 

a x 

We assume here a suitable discretization of the transverse plane with lattice points 
x. In Eq. ( |4.7|) , D is the total number of degrees of freedom, D = (N^ — ^)N X 
(i.e., the number of colors times the number of points in the transverse lattice). 
Recall that ovi, Xi-i are functionals of U^i so that Vi[a] depend on all the 
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variables «x, a 2 , ■ ■ -cti-i (c.f. Eq. ( |4.2| )). It is readily verified, by proceeding as 
indicated after Eq. ( |3.40| ), that the weight function ( [4.6[ ) is properly normalized: 



[da] W n [a] 



t d «] = nnn^)- 



(4.10) 



tax 



Since the inverse matrix \ is not immediately available, it is convenient to intro- 
duce an extra integration over auxiliary momenta irf(x) and transform Eq. (|4[ 
into 



W n [a] = / [d7r] expi-ie^ / a 

J i=1 Jx 

xexp{-6^ | / < X fi 



Ti ~ 1 / Ti^i-l 



(4.11) 



where [d7r] = J^ i=1 n a rice ( e /2 7r )d7r"(a;). It can then be easily verified by direct 
differentiation that W r [at] is solution of the following functional evolution equation 
which is the immediate translation of Eq. ( |3.48| ): 



d T W T [a] 



5 



[2 Sa b T (y) 



X±[a)W T [a)) -<£[a)W T [a] 



(4.12) 



As expected, this is the same as the renormalization group equation in the a- 
representation, Eq. Q2.16Q . 



4.2 The Fokker- Planck equation in the space of [/-fields 

Let us now define Z n [U] as the probability density for the random variable to take 
the value U at time r = ne, knowing that it is U at time 0. This is the analog 
of the probability density P(x,t) introduced in Sect. |3|. The normalization and 
initial conditions are written as 

J dfi(U)Z T [U] = 1 Z [U] = 5 G (U,U ), (4.13) 

where d[i(U) is the Haar measure and Sg(U, V) is the group-invariant 5-function 
(see Appendix). 

According to the discussion in Sect. 2.3, the probability Z n [U] can be calcu- 
lated from the weight function W n [a], via Eq. ( |2.19|) (for more transparency, we 
suppress space coordinates) : 

Z n [U] = J [da] W n [a] 5 G (U\ rf n [a]), (4.14) 

where Ul is given in terms of the a^'s by Eq. (|4.2p . (In fact, one could as well 
replace U\ by in the r.h.s. of this equation since lf n [a] ensures that a-,- = 
for any j > n). In continuous notations, Eq. (|4.14 ) becomes Eq. ( p.!9| ). 
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By exploiting the factorized structure of W n [a], Eq. (JO), one can verify that 



Z n [U] satisfies a recurrence formula similar to Eq. ( |3.26| ), that is, 

Z n [U] = jMUn-i) Z e [U\U n ^\ Z n _^J n _ x \ (4.15) 

and the elementary propagator Z e [U\U n -i] is given by an equation analogous to 
Eq. ( gig ): 



Z e [U i \U i . 1 ] = jldoc^VMSaiul^Ul^. (4.16) 

Note that for fixed C/j_i, the probability Vi[a] is in fact a function of aij given by 
Eqs. ( |4.7[ ) and ( [4.8[ ). We expect Z e [[/j|[Zj_i] to differ from unity by terms which 
vanishes as e — > 0. In fact the deviation is of order e, as we now show. To proceed, 
we note that the weight function Vi[a] is a Gaussian of width e in the variable 
ea, and we perform an expansion of the group (^-function in Eq. ( (4.16|) up to 
quadratic order in ec^ ~ sfe. To perform this expansion, we apply the identity 
( | A . 7| ) derived in Appendix and obtain, in direct generalization of Eq. ( |3.27| ), 



SGfuJ^Ul^Safe-^U}, U}_ x ) = e~^ a 5 G (ul U}_^ 

«(l - eay a + ^a^V a V b )s G (u}, U}_ x ) , (4.17) 

where the Lie derivatives act on £/j (not on Ui_ x ). By inserting this expression in 
Eq. ( |4.16|) and performing the integral over «j, we finally arrive at 



Z e [U\U n _ x ] = [l + e(-o a n _ x V a + \xtiV a V h )}8 G (u\ U\_ x ) + <D{e" 2 ). 

(4.18) 

The steps needed to obtain the equation satisfied by Z T [U] in the continuum 



limit are now identical to those leading to the Fokker-Planck equation (HT 
Reinstating coordinate dependence, one gets 



d T Z T [U] = -V a x V b y xty[U}Z T [U} - V a x a a x [U]Z T [U], (4.19) 



1 

2 

which is Eq. (|]|). 

On a more formal level, this equation could have been simply obtained by taking 
the time derivative of Eq. (|2.19|) relating Z T [U] and W T [a]. Noting that only W T 
depends on time in the r.h.s. of this equation, and using Eq. ( |4 . 1 2| ) for d T W T , one 
obtains: 

oo 

d T Z T [U] = J [da]d T W T 5 G (u\Te° driav U^ 

= / [da}W T [a] llx^la}— — ^— — - + a%[ a ]—L-\ 8 G (u\ Te'o u£) 

(4.20) 
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In the second line, we have performed integrations by parts w.r.t. a, and then re- 
stricted the integral in the exponent of the Wilson line to rj < r, taking advantage 
that the weight function VF T [a] imposes a" = for any r/ > r (c.f. Eq. ( |4.6| )). We 
then observe that the functional derivatives acting on this particular Wilson line 
have the same effect as the Lie derivative V a . Setting = Texp(i j^drj a v )Ul, 
we have indeed (see Eq. QA..6Q ): 



V\x)=it a V\y)5^{x-y)=V a y V\x). (4.21) 



5a^(y) ' ' y 

This allows us to trade the functional derivatives in Eq. (|4.20|) for Lie derivatives, 
e.g., 

A-S^UKTe^U') = V° v 6 a (U\V<) = -V^Te'^E/j). (4-22) 
and thus to immediately recover Eq. (|4.19|) . 



4.3 Path integrals in the space of ^/-fields 



We now return to Eq. ([4.14]) and use the known expression ([4.11|) for W n [a] to 
derive a path-integral representation for Z n [U] directly in terms of the C/j's, rather 
than the a^s. 

We proceed as usual, divide the rapidity interval r into n small elements of length 
e, and use repeatedly Eq. ( [4.15|) to deduce that 

/n— 1 
n^(^)] Z e [U\U n - X ] Z € [U n -x\U n - 2 ] ■ ■ ■ Z t [U x \U Q \. (4.23) 



Next, we define 



U\Ui-x = e 



ieuf = 2tr fhi(U}Ui-i). 



(4.24) 



The uiiS are uniquely determined as long as the matrix ufoi-i is sufficiently close 
to the unit matrix, which will be the case. 

Consider then the elementary propagator Z e [Ui\Ui-i] as given by Eq. ( |4.16| ). Using 
the parametrization ( |4.24j) , one can express the group 5-function in Eq. ( |4.16| ) as 
an ordinary ^-function for the random variable a" over which we need to integrate 
(see Eq. (|A.19|) in Appendix): 



\euj?t a ieaH c 



) = 8{euJi - eoti). 



(4.25) 



It is then straightforward to perform the integral over a" in Eq. ( |4.16|) , and obtain 
(detxi_i) 1/2 



i-1 



(2ne) D / 2 
x exp 



x 

(^-^Ux^C^-^t}, (4-26) 
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where it is understood that ujf(x) is a function of Ui(x) and Ui-i(x) obtained 
by inverting the relations ( |4.24| ). Eq. ( |4.26| ) is the analog of Eq. fl3.15| ), but it is 
actually more complicated than the latter because uji is non-linearly related to 
the neighboring group elements C/j and U{-\ (cf. Eq. ( |4.24|) ), as required by the 
geometry of the group manifold. 

Eq. ( [4.26| ) implies that, typically, euo a ~ ^/e, which is consistent with Eq. ( |4.25| ) 
and the previous estimates for «j (cf. Eq. ( |4.3|) ). This shows that, in order to 
control Eq. ( f4.26p to the accuracy of interest (i.e., up to and including terms of 
order e), it is enough to evaluate the measure factor fi(eui) to quadratic order in 
eoji. An explicit calculation to be presented in Appendix [A] yields (cf. Eq. ( |A.14[ )): 

W^iT 1 = 1 + e 2 ^ufuf + O^' 2 ) = e^« + O^' 2 ). (4.27) 



We are now in a position to write Z T [U] as a path integral: 

Z T [U] = I [dn(U)\ exp{-A[U}}, (4.28) 

where 

n-1 

[dfi(U)] = n {det Xl -iy 1/2 (2ne)- D / 2 [^(e^)]- 1 d^(^) (4.29) 



and 



n „ 



Again, one can trade the matrix \ 1 for % itself by introducing an additional 
integration over auxiliary momentum variables: 

/n— 1 „ n 

Y[dfi(U t ) / IJdTr* e~^ c ^ u \ (4.31) 
i=\ J i=i 



with 

U] = - e -§[ ^1 + i / < ( Wi - + ~ / < x£ i *f , (4-32) 

Jx Jx ^ Jx,y 

and we have also used Eq. ([4.27]) . To reintroduce the coupling constant g in the 
previous formulae, it is enough to insert a factor 1/g in the r.h.s. of Eq. (fOjj) 
for u}£. 

The role of the measure term (i.e., the first term in the r.h.s. of Eq. ( |4.32|) ) should 
be quite clear with the following remark: Assume that we decide to evaluate the 
path-integral ( |4.28| ) by using the "relative angles" u>i (in the sense of Eq. ( |4.24| )) 
as the integration variables. Then, as obvious from Eq. ( f4.29|) , the measure terms 
simply drop out: [^(tUj)]" 1 d/i(£7j) = dcjj, and Eq. ( |4.28| ) reduces to Eq. ( [1.14| ) 
with W n [a] of Eq. 
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5 Conclusions 



We have established the relationships between the two forms of the renormal- 
ization group equation towards small x which are available in the literature. 
These have been obtained in different approaches, and, in spite of some for- 
mal similitudes, they rely on different mathematical structures: the U- and the 
a-representations of this paper. What we have shown here is that these two 
representations are equivalent and describe the same basic statistical process: a 
random walk in the space of Wilson lines, with In 1/x playing the role of "time" , 
a that of the "velocity" , and U being the "position" on the group manifold. 

The Langevin description of this random walk provides a natural bridge between 
the two representations: At step i, the gauge field «j is randomly chosen according 
to the gaussian probability distribution (}4.3|) , then used to update the SU(N) 
group element Ui, cf. eq. (p~T|), which is in turn used to compute the probability 
law for selecting the "velocity" a i+ i at step i + 1, and so on. This Langevin 
process provides an algorithm for the quantum evolution and stays at the basis 
of a numerical study of the renormalization group equations which is currently 
in progress K| . Note that the results of such Langevin simulations will not only 
give access to the x dependence of cross sections like the j*A cross sections of 
Eq. ( |1.2|) , they also provide the basis for an extension to the determination of 
initial conditions in AA collisions with a calculated x dependence along the lines 
of §,§. 



In addition to its suitability for practical implementations, the Langevin formu- 
lation brings also important conceptual clarifications, as it most prominently ex- 
hibits the statistical nature of the underlying physics: The renormalization group 
at small x is driven by successively integrating out gluonic fluctuations within a 
small slice in In 1/x around a static (and typically large) background field, which 
represents the quantum modes integrated in previous steps, and which are ef- 
fectively "frozen". The ensuing coarse graining in the longitudinal direction is 
accompanied by a loss of quantum coherence, allowing eventually a description 
in terms of a classical stochastic process. 

The particular structure of the stochastic equations which emerge from the renor- 
malization group has allowed us to avoid discretization ambiguities which gen- 
erally arise for random walks in inhomogeneous media. We have utilized the 
so-called Ito discretization, which allows one to construct the new variables in 
a "causal" manner from the knowledge of the old ones, and we have verified 
that, with this choice, the Langevin description reproduces indeed the expected 
renormalization group equations — which here emerge as Fokker-Planck equa- 
tions — in the continuum limit. Once things are properly formulated within 
one discretization, other prescriptions can be used as well, and we have briefly 
elaborated on some alternatives. 

The discretization that we have used has also the virtue to greatly facilitate 
the construction of path-integral representations for the solutions to the Fokker- 
Planck equations. Such path-integrals can be the starting point for alternative 
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numerical implementations. It would be interesting in this context to explore the 
possibility to rewrite this path-integral in a form which is local in the transverse 



coordinates, possibly after introducing auxiliary fields, as in Ref. [29]. But path- 



integrals can be also useful for analytical studies, and could give us more insight 
into the qualitative features of the non-linear physics at small x. Although to 
some extent formal, the explicit path-integral solutions that we have obtained 
have the merit to exhibit important structural properties, like non-locality (in 
both longitudinal and transverse coordinates), non-linearity, and colour struc- 
ture. Less formal analytic results could be obtained by further evaluating the 
path integrals, e.g., via saddle-point or mean- field approximations. In particular, 
the self-consistent mean-field approximations discussed in Ref. || can be easily 
reformulated within the present approach. 
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A Calculations on Lie group: useful formulae 

We collect in this appendix a few identities which are useful in calculations on 
group manifolds (see, e.g., Refs. ^6|,01 f° r more details). 

We consider SU(N) matrices U. These depend on N 2 — 1 real parameters u a , and 
can be written in the exponential form as: 

U(u) = e~ iujata = e-' lUJ U\uj) = e iUJ (A.l) 

where the t a 's are the infinitesimal generators normalized as Tr(t a t b ) = dR5 ab , 
with d,R the dimension of the representation considered. When needed, we shall 
distinguish the adjoint representation by a tilde: U(u) = e~ luJ . 

Consider first the case where these parameters depend on a single variable, the 
"time" t. When r varies, the matrix U T describes a curve on the group manifold. 
Consider a variation of the matrix U when r increases by a small amount dr. We 
can write: U T+ d T = U T + dU T , and note that U^dU T = U^.U T+ d T — 1 is an element 
of the Lie algebra that we can write as UJj.dU T = —iT a a a dT where a a (r) are the 
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components of a vector tangent to the curve. Clearly dU T = —iU T T a a a dr, or 
equivalent ly 

dU j T = ia a T a Uldr (A. 2) 

which allows us to write £/J as a time ordered exponential: 

U\ = Texp{* f dr]a a (rj)T a }Ul (A.3) 
Jo 

Note that it is the matrix W, rather than U itself, which enters our discussion in 
the main text. Note also that a a ^ du a /dr. 

Consider next an arbitrary function of C/J. We may write: 

^1 = (»T^) y . * /[C/t] = a a V<7(^), (AA) 

t/y 

and therefore 

/(C/t) = Texp y\ V a a ( V ) V a }} (A.5) 

The operators V a introduced above, commonly referred to as Lie derivatives, 
satisfy the following relations: 

V a [/ = -iU T a , V a U ] = iT a U ] , [V a , V b ] = f abc V c , (A.6) 

and form therefore a representation of the Lie algebra. 

If one chooses a to be constant in eqs. ( |A.3|) and ( |A.5|) , one gets the following 



useful formula: 

fie^Ut) = e^ a f(Ul) (A.7) 

On the group manifold, one may define a metric by ds 2 = Tr (dUdW) or, in 
terms of the group parameters u a : 

^ = g ^ = Tr^^ = (LL%, (A.8) 

where we have set 

= -iL ab {uo)T\ (A.9) 

and L* fe = L^. From this one gets easily the Haar invariant measure 

d/i(U) = fi{uj)U a duj a n(u) = y/\detg{u)\. (A.10) 



27 



To calculate L ab we use: 



Note the emergence of adjoint-representation matrices Co in the r.h.s.'s of the 
above formulae. Since matrices in the adjoint representation are real, {e~ lsw ) a b = 
(e lsu} )ba)- One can then write 

(L(u)L\u)) ab = f Q dsj\\ (e^) ab = Q J^dse^ 2 , (A.12) 



and obtain ix{ui) = | det(j7" a fe(cj))| with 

= 1 f ds lei-*) - v { - )n ^» - (^m\ , A 131 



This expression allows in particular the computation of to quadratic order 
in w, as required in Sect. 4.3: 

= det J(lo) = e Trln ^ } = e~> a " a + 0(lo 3 ). (A.14) 
Integrals over the group manifold are typically of the form: 

J dn(U)f(U) = j U a dcu af x(co)f(e-^). (A.15) 
A particularly useful tool is the group 5-function defined by 

MU)f[U]5 G [U,V}=f[V] (A.16) 



with obvious invariance properties, like 5 G [U, V] = 5 G [UV\l] = 8 G [1,WV], 
which follow directly from the invariance of the Haar measure dfj,(U). It will 
useful to relate this group 5-function to the ordinary 5-function. To do so, con- 
sider the integral: 

f(e- i0 ) = jdpt(U)f(U)S G (U : e- i0 ) 

= Ju a doj af i(u;)f(e- i<jJ )5 G (e- iUJ ,e- i0 ) (A.17) 

On the other hand, we have 

f(e~ id ) = J n a dcu a /(e- iW ) 5{u - 9), (A.18) 

with 5(u — 9) = U a 5(uj a — 9 a ). Combining these two equations, we obtain the 
following identity: 

/i(u)5 G (e- i0 , e~ iUJ ) = 5{uj - 9) (A.19) 
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Further useful formulae involve Lie derivatives. First we may use e.q. flA.llj) 
together with the properties ( |A.6| ) of Lie derivativs to obtain: 



dU 



L ab V b U 



L ab V b Ul 



(A.20) 



From these relations, and again e.q. ( A.ll|) , one easily deduces the following useful 
identities: 







V a f(U)\ u = 1 = ^- a f(U(u J ))\ ul=0 



(A.21) 



\ (V"V b + V fc V") f(U)\ u=1 = /(C/(o;))L =0 . (A.22) 

By using these fromulae and the relation obtained above for the group ^-function, 
one can also derive the following identities: 

»(u)V a 6 G {U(u), 1) = ^S(ou) , (A.23) 



^) x ab V a V b 5 G (U(uj), 1) = x^^-^Sico) , (A.24) 
for any symmetric x ab ■ 

All the previous formulae are easily extended to the case where the matrices U 
are fields U(x) in the transverse coordinates x. Thus for instance, Lie derivatives 
become 

Kfl U]=(-iU { x)t% w l ( - j flUl (A.25) 

We have also: 

V a x U y = -iU y t a 5 xy , V a x U y = it a U y 5 xy , [V^, V* ] = f abc V c y 6 xy , (A.26) 
where 5 xy = 5 (2) (x - y). 
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